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abstract 

Soliton motion in some external potentials is studied using the nonlinear Schrodinger 
equation. Solitons are scattered by a potential wall. Solitons propagate almost 
freely or are trapped in a periodic potential. The critical kinetic energy for 
reflection and trapping is evaluated approximately with a variational method. 
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The nonlinear Schrodinger equation has been intensively studied as a soli- 
ton equation, which is derived as an envelope equation by several perturbation 
methods pQ. Solitons in optical fibers are typical ones described by the non- 
linear Schrodinger equation 0. Recent successful observation of solitons in 
Bose-Einstein condensates (BECs) raises the theoretical consideration of this 
equation again 3 . The equation is called the Gross-Pitaevskii equation in 
the research field of BECs. In the BECs, matter waves obey the nonlinear 
Schrodinger equation with external potentials. For example, external magnetic 
fields are applied to confine the BECs. Interference patterns of laser beams 
generate external periodic potentials. Solitons in external potentials such as 
harmonic potential, periodic potential and stepwise potential have been studied 
by several authors 010311)] ■ But solitons in a potential wall or potential well were 
not studied in detail. In this letter, we will study the scattering and trapping of 
solitons by a potential wall (or well) and a periodic potential. In particular, we 
focus on soliton behaviors like classical mechanical particles, and evaluate the 
critical energies for the transmission and the trapping. 

The nonlinear Schrodinger equation with a potential term is written as 

where 4>(x,t) is the wave function, U(x) is an external potential. Without the 
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nonlinear term, this equation is equivalent to the Schrodinger equation with 
mass m = 1 and % = 1. When U(x) = 0, the nonlinear Schrodinger equation 
has a soliton solution 

<j>(x, t) = Asech{A(x - vt)}^*-^-™* , (2) 

where v = p and lu = -A 2 /2 + v 2 /2. 

Firstly, we investigate the nonlinear Schrodinger equation in a simple rect- 
angular potential, that is, 

U(x) = Uo, for x < x < x + d, U(x) = 0, for x < x n , x > x + d. (3) 

If Uo is positive (negative), U(x) represents a potential wall (well). We consider 
first the case of Uo > 0. If the nonlinear term is absent in cq. (1), this is a typical 
potential problem of the quantum mechanics. The transmission coefficient T of 
matter waves is expressed as 



T = 



1 + 
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for E > U , T = 
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for E < Uo, 
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where E denotes the kinetic energy E — v 2 /2, and a and a' are defined as 
a = y/2jE - Uo) (for E > Uo) and a' = ^2(Uo - E) (for E < Uo). We study 
numerically the scattering of a soliton by the potential wall. Figure 1(a) displays 
the time evolution of a soliton passing through a potential wall with Uo — 0.5 
and d=l. The system size is L = 100 and the potential wall locates at x = x = 
L/2 = 50. The initial velocity of the soliton is v = 0.75 and the amplitude is A = 
1. The soliton is split into two pulses by the potential wall. The transmission 
coefficient is 0.527, which is calculated as T = J L ^ 2 \(f>\ 2 dx/ J Q \<p\ 2 dx after the 
scattering. Figure 1(b) displays the transmission coefficient as a function of the 
initial velocity v. The dashed curve is the transmission coefficient eq. (4). The 
transmission coefficient of a soliton is close to the coefficient eq. (4) for v > 0.9, 
which might imply that the effect of the nonlinear term becomes relatively 
weaker as the kinetic energy is increased. The transmission coefficient decreases 
more sharply than the quantum mechanical particle as v is decreased. 

The transmission coefficients as shown in Fig. 1(b) change more and more 
sharply as Uo is decreased. As Uo becomes smaller, such a splitting or ad- 
ditional radiations are hardly observed and the soliton keeps its form and 
behaves like a particle. The coefficient jumps almost from to 1 at a cer- 
tain critical velocity for Uo < 0.05. To understand the transition between 
the wavy behaviors including the splitting and the particle-like behavior, we 
have performed numerical simulations for another initial condition <j>{x, 0) = 
Asech{A(a; — L/2)}, that is, the soliton is initially set at L/2 and the ini- 
tial velocity is 0. The potential wall is located in [L/2 — d/2.L/2 + d/2}. 
The mirror-symmetry around L/2 is maintained during the time evolution. 
For large Uo, the soliton is split into two pulses as shown in Fig. 2(a) for 
Uq = 0.67, A = 2, d = 1 and L = 80 owing to the repulsive potential. How- 
ever, the soliton cannot be split for small Uo owing to the attractive interaction 
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Figure 1: (a) Time evolution of |0(a;,t)| for Uo = 0.5,^4 = 1 and v = 0.75. 
A soliton is splitted into two pulses by the potential wall of d = 1 located at 
x = L/2 = 50. (b) Transmission coefficient T as a function of the initial velocity 
v. The dashed curve is the transmission coefficient for the quantum mechanics 
expressed by eq. (4) 

by the nonlinear term. The soliton keeps its one pulse structure as shown in 
Fig. 2(b) for U = 0.66, A = 2, d = 1 and L = 80. The critical value of U 
depends on the initial amplitude A as shown in Fig. 2(c). Equation (1) has at 
least two invariants, that is, the total mass N = J \4>\ 2 dx and the total energy 
E = jf{l/2\(t) x \ 2 - 1/2|</>| 4 + U(x)\4>\ 2 }dx. If the soliton is split into two coun- 
terpropagating solitons: (f>(x, t) = A'sec\i{A'(x - vt - L/2)}e l ^ x - vt - L / 2 ^- lu}t + 
A'sech{A'(x + vt- L/2)}e- tp ( x+vt - L / 2 ')-^ t , the amplitude A' and the momen- 
tum p satisfy A' = A/2 and 2p 2 A' = -1/4A 3 + 2U Ata,nh(Ad/2) owing to 
the conservation laws. If Uo < A 2 / {8 tanh(Af/2)}, there is no real solution 
for p, that is, the splitting cannot occur and the one-pulse structure is main- 
tained. Above the line Uo — A 2 / {8 tanh( Ad / 2)} , the potential energy is used to 
split a soliton against the attractive interaction. A soliton behaves like a par- 
ticle even for more general initial conditions, if Uo is sufficiently smaller than 
A 2 /{8tanh(Ad/2)}. 

We focus on the particle-like behavior below. We return to the numerical 
simulations of the first type. Figures 3(a) and (b) display the time evolutions of 
a soliton passing through a potential wall of A = 1,Uq = 0.04 and d = 1. The 
system size is 50 and the potential wall locates at xo = L/2 = 25. The soliton 
is almost completely reflected when the initial velocity is v = 0.18, and almost 
completely transmitted for v = 0.2. The critical velocity for the soliton to pass 
through the potential wall is 0.1925, and the critical kinetic energy is 0.0185, 
which is smaller than the critical value E c = Uo — 0.04 for a particle obeying 
the classical mechanics. Figure 3(c) displays the critical kinetic energy as a 
function of the amplitude A of the soliton. As the amplitude A is increased, the 
critical energy approaches Uo- A soliton with smaller amplitude can pass easily 
through the potential wall. This may be interpreted as coherent tunneling. 
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Figure 2: (a) Time evolution of \(j)(x,t)\ for U = 0.67, A = 2 and L = 80. 
A soliton is splitted into two pulses by the potential wall of d = 1. (b) Time 
evolution of |(/>(a;,i)| for Uq — 0.66,^4 = 2 and L — 80. A soliton keeps its 
one-pulse structure, although it is breathing, (c) Critical value of Uo for the 
splitting of the soliton as a function of the initial amplitude A. The dashed line 
is the theoretical one Uq — A 2 /{8 t&nh( Ad / 2)} 



The nonlinear Schrodinger equation is rewritten with the Lagrangian form: 
d/dt(5L/54> t ) = 5L/S(f>, where L = f™ oo dx{i/2(<t> t <j>* - (j>*(j>) ~ l/2\cj> x \ 2 + 
l/2\<j)\ 4 — U(x)\<j>\ 2 }. If a solution to the nonlinear Schrodinger equation is 
approximated by 

<j>(x,t) = A{t)sec\y{{x - C(t))/a(t)}e jp(t)(K -« (t)) - Mt , (5) 

the time evolution of a(i),£(i) and p(t) can be derived from the effective La- 
grangian L e f f = Ldx using 

L eff = 2A 2 aptt + 2/3A 4 a - l/3A 2 /a - A 2 ay 2 - 2A 2 aU eff , 

where U e ff — ?7o/2[tanh{(a;o+(i— ^) /a}— tanh{(xo— ^)/a}]. Equations of motion 
for the variational parameters qt = £,p and a are given by d/ dt(dL e f f / dq\) — 
dL e ff/dqi under the condition N = \cj)\ 2 dx = 2A 2 a = const. The results 
are 

d£_ dp _ dU eff da =( (( ,, 

dt Pi dt ds, ' dt { ' 

and A = 1/a = N/2. The soliton obeys the Newton equation with an effective 
potential U e ff. The peak amplitude of the effective potential is Uq tanh( Ad / 2) . 
The critical kinetic energy is therefore expected to be E c = Uq tanh( Ad / 2) , 
which was drawn with the dashed curve in Fig. 3(c). Good agreement of nu- 
merical results with this approximation is seen in Fig. 3(c). 

Secondary, we consider the case of negative Uq. The potential acts as a 
potential well. We study the scattering of a soliton by the potential well. Com- 
plicated scattering including splitting is observed when \Uq\ is large, as is the 
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Figure 3: (a) Time evolution of \</>{x,t)\ for U = 0.04, A = 1 and v = 0.18. (b) 
Time evolution of \<j>(x, t)\ for Uo — 0.04, A = 1 and i> = 0.2. (c) Critical kinetic 
energy in the potential wall of Uq = 0.04 as a function of A. The dashed curve 
is a theoretical approximation 0.04 tanh(^4/2) 



case of positive Uq. We focus on the particle- like behavior, which is observed 
when \U \ is sufficiently small. Figure 4 displays the time evolutions of solitons 
for a potential well with Uq = —0.04 and width d = 1. The initial amplitude 
is A = 1 and the initial velocities are v — 0.035 and 0.015. For relatively large 
v, the soliton passes over the potential well like a classical mechanical particle 
obeying eq. (6). The velocity increases when the particle passes over the poten- 
tial well, but its initial value is recovered after the passage. However, the soliton 
is trapped by the potential well for sufficiently small v as shown in Fig. 4(b). The 
critical value of v for the trapping is approximately v c ~ 0.0225. In Fig. 4(a), 
a soliton passes through the potential well, but the velocity is decreased from 
0.035 to 0.026 after the passage of the potential well. Some dissipative processes 
may occur at the potential well. Some small radiations may be generated and a 
part of the kinetic energy is changed into radiation energy. Another possibility 
is that some internal freedoms of a soliton such as breathing motion may be ex- 
cited and the kinetic energy decreases. To investigate such motion, we assume 
a solution of the form 

(j){x,t) = A{t)scch{(x _^( t ))/ a ( t )} e <P(t)(»-€(*))+to(*)logcoBh{(x-€(t))/a(t)}-ia,t i 

(7) 

where A, a,p, £ and a change in time, but N — 2A 2 a is conserved in time. The 
term ia(t) logcosh(.x — £(t))/a(t) is a kind of chirp term and this form of chirp 
term appears in a solution to the complex Ginzburg-Landau equation where 
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Figure 4: Time evolutions of |(/>(x, t)| in a potential well of Uq = —0.04 and 
d=l. The initial velocity is 0.035 (a) and 0.015 (b). 



some dissipation and filtering effects arc involved. Owing to the inclusion of 
this term, the width of the soliton changes in time and breathing motion be- 
comes possible. The choice of this form of chirp term is not a unique one. We 
have studied another form of chirp term ia(x — £(t)) 2 , but better agreement 
with direct numerical simulation was seen for the form eq. (7). The substitu- 
tion of eq. (7) to the Lagrangian and its spatial integration yield the effective 
Lagrangian 

L eff = N(p^ t -a t (l-log2)+aa t /(2a)-l/(6a 2 )-p 2 /2-a 2 /(6a 2 )+N/(6a)-U eff ), 

where U e ff — J7 /2[tanh{(x +d — £)/a} — tanh{(a;o — £)/°}]- The equations of 
motions for p(t), a(t) and a(t) are 

d£ dp _ dU eff 

dt ~ P ' dt~ d£ ' 

± ^ d(T -n n f 1 + (j2 N 9U 'fA 

dt ~ 3a' dt V 3a 3 6a 2 da J ' { ' 

The energy p 2 /2+U e f / is not conserved in this equation, and the width a changes 
in time. If the energy decays, the soliton may be trapped in the potential 
well. Figure 5 displays the time evolutions of by eq. (8) for v = 0.035 
and v = 0.015 with the dashed curves, and they are compared with the time 
evolutions of the peak position of the solitons shown in Figs. 4(a) and (b). The 
velocity of a particle changes from the initial value of 0.035 to 0.029 in the time 
evolution of eq. (8) in Fig. 5(a). The particle is trapped in the potential well 
and is oscillating in a region of 23 < x < 28 for the initial velocity v = 0.015 in 
Fig. 5(b). The trapping occurs for v < 0.019 in eq. (8). The dissipative effect 
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Figure 5: Comparison of the time evolutions of the peak position X p (t) of a 
soliton (solid line) and by Eq. (8) (dashed line) for v = 0.035 (a) and 0.015 
(b). 

seems to be slightly weaker in eq. (8) than in the nonlinear Schrodinger equation 
eq. (1). 

As the third case, we study soliton motion in a spatially periodic potential 
U(x) = — Uq cos(27rx). Initially, a soliton locates at the position of the potential 
minimum x = L/2 = 25, and the amplitude is A = 2. Figure 6 displays the time 
evolutions of solitons in the periodic potential with Uo = 0.2 for two different 
initial velocities v. A soliton is almost steadily propagating for v = 0.33 as shown 
in Fig. 6(a). For v = 0.28, a soliton propagates in an initial stage, however, it is 
trapped in the periodic potential as shown in Fig. 6(b). The critical velocity for a 
soliton to be trapped for A = 2 is approximately 0.312. Figure 7(a) displays the 
critical kinetic energy E c = as a. function of A. The solid curve represents 
the critical kinetic energy, and the rhombi represent critical points below which 
a soliton cannot go over the nearest potential hump. (For the parameter of v 
shown in Fig. 6(b), the soliton goes over several potential humps and then it is 
trapped.) The effective potential U e ff for the periodic potential is expressed as 

f°° , , ,9 , -C/ cos(27rf)7r 2 a 

U eff = {-U cos(2irx)}\(l)\ 2 dx/N = 



sinh(7r 2 a) 

A particle, which obeys eq. (6) with the effective potential t/ e // and starts from 
a position of potential minimum, exhibits a transition from a trapped state to a 
propagating state. The critical initial kinetic energy is 2[/o7r 2 /(^4sinh(7r 2 /^4)), 
which is drawn by the dotted curve in Fig. 7(a). The critical curve is a good 
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Figure 6: Time evolutions of \<p(x, t)\ in a periodic potential U(x) — 
— 0.2 cos(27rx). The initial amplitude A is 2, and the initial velocity is v = 0.33 
(a) and 0.28 (b). 



approximation for the critical values marked with rhombi. This is probably 
because eq. (6) is a good approximation to the short-time evolution of eq. (1), 
where dissipative effects are negligible, and which determines whether the soliton 
can pass over the first neighboring hump. In the region 1.5 < A < 3, the 
critical curve 2Uq-k 2 /(Asinh(7r 2 /A)) is rather different from the true critical 
curve denoted by solid line. If the width of a soliton 2/ A is comparable to 
the wavelength of the periodic potential, the interaction between the soliton 
and the potential will be enhanced, which may be an origin of the deviation. 
The trapping phenomena accompanying the dissipation of kinetic energy can 
be expressed with eq. (8) with U e ff = — Uq cos(27r£)7r 2 a/ sinh(7r 2 a). The time 
evolution of by eq. (8) for A = 2 and v — 0.2785 is displayed in Fig. 7(b) with 
the dashed curve. The solid curve in Fig. 7(b) is the time evolution of the peak 
position X p (t) of the soliton for A = 2 and v = 0.28. Similar time evolutions 
including damping of kinetic energy, trapping and chaotic motions are observed 
for £(t) and X p (t), although the initial velocities are slightly different. Since the 
number of degrees of freedom in eq. (8) is four, it is not surprising that chaotic 
motion appears in eq. (8). The dashed curve in Fig. 7(a) is the critical line, 
above which steadily propagating solutions appear in eq. (8). There is fairly 
good agreement between the solid line and the dashed line, although the dashed 
line locates slightly below the sold line in the region 1.5 < A < 3. 

In summary, we have performed numerical simulations for the scattering and 
the trapping of nonlinear Schrodingcr solitons in some external potentials. The 
transition from the reflection to the transmission by a potential wall depends 
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Figure 7: (a) Critical kinetic energy for the propagation in a periodic potential 
U(x) = — 0.2 cos(2-7ra) as a function of A. The solid line denotes the critical 
kinetic energy, above which a soliton propagates almost steadily. The rhombi 
denote critical values, below which a soliton cannot go over the nearest potential 
hump. The dotted curve is a theoretical line OAw 2 / (As'mh(ir 2 /A)). The dashed 
line denotes a critical kinetic energy above which a particle obeying eq. (8) 
propagates almost steadily, (b) Time evolutions of the peak position X p (t) by 
eq. (1) (solid line) and by eq. (8) (dashed line). 

strongly on the soliton amplitude. We have also found the trapping of solitons in 
a potential well and a periodic potential. The critical energy for the transitions 
can be approximately evaluated using eq. (6) or (8) derived from a variational 
method. 
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